arXiv:1503.03395vl [astro-ph.SR] 11 Mar 2015 


Astronomy & Astrophysics manuscript no. radial3 arxiv © ESO 2015 

March 12, 2015 


Stopping Frequency of Type III Solar Radio Bursts in Expanding 

Magnetic Flux Tubes 

Hamish A. S. Reid and Eduard R Kontar 
SUPA School of Physics and Astronomy, University of Glasgow, G12 8QQ, UK 


ABSTRACT 

Aims. Understanding the properties of type III radio bursts in the solar corona and interplanetary space is one of the best ways to 
remotely deduce the characteristics of solar accelerated electron beams and the solar wind plasma. One feature of all type III bursts is 
the lowest frequency they reach (or stopping frequency). This feature reflects the distance from the Sun that an electron beam can drive 
the observable plasma emission mechanism. The stopping frequency has never been systematically studied before from a theoretical 
perspective. 

Methods. Using numerical kinetic simulations, we explore the different parameters that dictate how far an electron beam can travel 
before it stops inducing a significant level of Langmuir waves, responsible for plasma radio emission. We use the quasilinear approach 
to model self-consistently the resonant interaction between electrons and Langmuir waves in inhomogeneous plasma, and take into 
consideration the expansion of the guiding magnetic flux tube and the turbulent density of the interplanetary medium. 

Results. We find that the rate of radial expansion has a significant effect on the distance an electron beam travels before enhanced 
leves of Langmuir waves, and hence radio waves, cease. Radial expansion of the guiding magnetic flux tube rarefies the electron 
stream to the extent that the density of non-thermal electrons is too low to drive Langmuir wave production. The initial conditions of 
the electron beam have a significant effect, where decreasing the beam density or increasing the spectral index of injected electrons 
would cause higher type III stopping frequencies. We also demonstrate how the intensity of large-scale density fluctuations increases 
the highest frequency that Langmuir waves can be driven by the beam and how the magnetic field geometry can be the cause of type 
III bursts only observed at high coronal frequencies. 
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1. Introduction 


Type III r adio bursts observed in the meter range were initially 
classified ( [Wild & McCready||19501 ) as an impulsive radio fre¬ 
quency signal that drifted from high to low frequencies with 
time. They were given the name ‘type III’ because the rate they 
drift in frequency with time (or drift rate) was higher than ‘type 
I’ and ‘type II’ bursts. Type III bursts are one of the most fre¬ 
quently observed impulsive electromagnetic signals driven by 
the Sun and are produced by accelerated electrons beams trav¬ 
elling at near-relativistic energies through the plasma of the 


Dulk 


solar system (see e.g. |Goldman| 1983|JSuzuki & Dulk||T 985 ; 


2014 


[M5l IMelrosell 19901 |Muschiitti|| 19901 |Reid & Ratcliffel 


for reviews). Understanding which properties of acceler¬ 
ated electron beams and the solar wind plasma give rise to cer¬ 
tain features in type Ills allows us to use these transient burst 
properties as a diagnostic tool. 


The escaping radio emission responsible for type III bursts 
can be broadly described as two stage process (first suggested 
by [Ginzburg & Zhelezniakov||1958| ). The first stage is a two- 
stream instability between the electron beam and the background 
plasma that induces a high level of Langmuir waves. The second 
stage is wave-wave interactions between Langmuir waves and 
either ion-sound waves for fundamental emission or oppositely 
propagating Langmuir waves for harmonic emission. The cur¬ 


rent theory has been refined by many authors (e.g. Zheleznyakov 


& Zaitsev 1970 

Smith! 

1970; Zaitsev et al. 1972f Smith et al. 

1976 

Melrose 

1980b; 

Goldman 

1983; Dulk 1985; Melrose 

1987 

Kontar et al. 1998; Mel’Ni 

k|[2003j Ledenev et al.p004) 


and continues to develop with the advent of numerical simu- 


lations (e.g Takakura & S 

libahashi 1976; Magelssen & Smith 

1977 

, Grognard 1985| Ro 

Dinson 1992; Kontar 2001d; Li et al. 

2008 

, Kontar & Reid 2009 

Tsiklauri 2010 ,2011; Reid & Kontar 

2013 

, Ratcliffe et al.|2012[ 

Li & Cairns 12014; Ratcliffe & Kontar 

2014 

). 


Type III bursts cover a wide range of frequencies. The start¬ 
ing frequency of type III bursts varies from hundreds of MHz, 
and rarely GHz, down to 1 MHz and below. The stopping fre¬ 
quency of type III bursts varies over a even wider frequency 
range from the 100’s to 10’s MHz for type III bursts that origi¬ 
nate only in the corona, all the way down to 10’s kHz for type III 
bursts from interplanetary space. The properties of an electron 
beam and the background heliospheric plasma govern at what 
frequencies the type III radio emission starts and stops being 
observed. Whilst recent studies have investigated what is im¬ 
portant in determining the starting frequency of type III bursts 
( Reid et aL||2011| |2014| ), the processes affecting the stopping 
frequency of type III bursts has never been systematically in¬ 
vestigated. Figure [T] shows an example of many type III bursts 
observed before and after a large flare. Although electron beams 
associated with the flare reached frequencies around 0.04 MHz, 
the majority of other type III bursts stopped at higher frequen¬ 
cies. What plasma parameters dictate at what frequency a type 
III bursts stops is unclear. 

Previous observational studies on type III burst stopping fre¬ 
quencies by |Leblanc et al.| ( |1995| [1996) used radio data from 
the WAVES instrument onboard the WIND spacecraft and th e 
URAP experiment on the Ulysses spacecraft (Stone et al. 1992). 
They analysed type III bursts detected only by Ulysses between 


1990 and 1994 ( [Leblanc et al.||1995] ) and by both WIND and 
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16 FEB 2014 



tion of radio emis sion (previously suggested by |de G enouillac 
& Escande ( [1981 }). In the third, the density fluctuations sup¬ 
press the level of Langmuir waves excited by the electron beam 


(e.g. Nishikawa & Ryutov|19761 ) and thus cease the production 


of radio emission. That large scale density fluctuations suppress 
Langmuir wave growth has been discussed at length by many 


authors (e.g. Smith & Sime 

1979; Muschietti et al. 1985 Kontar 

200 lH Reid & Kontar 2010 

nUet al. 2011; Ratcliffe et al. 2012). 


Dulk et alT| ( [1996[ ) tested the hypothesis that radio wave pro¬ 


rogation governs the lowest type III frequencies and that the cut¬ 
offs are not intrinsic to the radiation mechanism. Using simulta¬ 
neous WIND and Ulysses (which was behind the Sun) observa¬ 
tions, they found a reasonable correlation coefficient of 0.51 but 
with a high degree of scatter. 28% of bursts they observed had 
the same stopping frequency to within a factor of 1.12, whilst 
50% of the bursts had the same stopping frequency to within a 
factor of 1.4. At the extremes, the stopping frequency of type III 
bursts differed by a factor of 5. The difference in cutoff frequen¬ 
cies suggests that some directivity or rather the blocking of radio 
emission between source and observer by density fluctuations 
can play a role in determining stopping frequency. However, it 
is not the dominant process because the stopping frequency of 
type III bursts can vary over many orders of magnitude in fre¬ 
quency. Lor example in the study by [Dulk et aL| ( [1996} ) they 
recorded stopping frequencies from 300 kHz down to 20 kHz. 
The large spread suggests that the plasma processes involved 
with the production of type III bursts largely determine the stop¬ 
ping frequency. 


Fig. 1 . An example of many type III bursts on the 16th Lebruary 
2014, detected by the WAVES instrument on the WIND space¬ 
craft at 14 - 0.02 MHz (Bougeret et al. 1995[) , the Nan 9 ay 


A more recent study using the WAVES instrument ( Bougeret 
et al.|2008]) onboard the STEREO spacecraft was carried out by 


Decametre Array at 80 -14 MHz (Lecacheux]2000) and Phoenix 
4 at 900 - 200 MHz ( [Benz et al.|2009| . The large type III burst 
at 09:30 was associated with a GOES Ml flare. We can observe 
type III bursts stopping at a variety of different frequencies. The 
periodicity observed below 1 MHz is an artefact. 


Krupar et al. ( j2014| ) on 154 type III bursts occurred between 


2007 and 2013. They found that the stopping frequency in 65% 
to 75% of events was at 125 kHz, the lowest frequency chan¬ 
nel measured. It is also important to mention that the measured 
stopping frequency is dependent upon the sensitivity of the radio 
receiver and/or the level of ambient signal from the background 
plasma and the galactic background. 


Ulysses detected below 1 MHz between 1994 and 1995 ( |Leblanc| 
|et al.|1996| ). Each type III burst was given a stopping frequency, 
estimated from the lowest observable frequency on the type III 
dynamic spectra. Repeating the process a number of times for 
the same type III bursts they estimated an error on their results 
of 12%. Typical stopping frequencies for type III bursts were 
observed between 100-200 kHz for weak type III bursts and 20- 
50 kHz for strong type III bursts. Neither study observed type 
III bursts below 9 kHz despite the plasma frequency at Ulysses 
being as low as 3 kHz during observations. They found a con¬ 
nection between the intensity of the type III burst emission and 
the stopping frequency. The larger the intensity of the type III 
burst, the lower the stopping frequency tended to be. 

There were three postulations by [Leblanc et al.| ( 1995| ) for 
what influences the burst stopping frequency. In the first, the 
electron beam is diluted due to expansion of the guiding mag¬ 
netic flux tube. The diluted electron beam becomes ineffective 
in forming the required magnitude of positive slope in velocity 
space required for generating Langmuir waves. In the second, 
the level of density fluctuations or ion-acoustic waves required to 
convert Langmuir waves into electromagnetic waves decreases 
as a function of distance from the Sun and throttles the produc- 


The necessary condition for a type III radio burst is a pop¬ 
ulation of coronal accelerated electrons that are able to drive 
Langmuir waves with access to open magnetic field lines. If a 
population of accelerated electrons is able to drive Langmuir 
wave generation, what causes this Langmuir wave generation to 
cease and how does it vary with distance from the Sun? To an¬ 
swer these questions, we have simulated the outward propaga¬ 
tion of a solar electron beam injected in the low corona. We self- 
consistently modelled the wave-particle interactions that govern 
the production of Langmuir waves through the beam-plasma in¬ 
stability. We also modelled collisional effects, wave refraction in 
an inhomogeneous background plasma, the expanding magnetic 
field lines of the solar corona and solar wind. The model and the 
initial condition of the injected electron beam are all described 
in Section [2] Using Langmuir wave energy density as a proxy 
for the level of radio emission generated by an electron beam, 
we considered a number of different effects that govern when 
the production of Langmuir waves will cease: the expansion of 
guiding magnetic flux tube in Section[3] the level of density fluc¬ 
tuations in Section |4| and the initial beam parameters in Section 
[5] Our results are discussed in Section [6] 
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2. Model equations, initial conditions and 
approximations 

2.1. Model 


Each dimension of the source function is characterised by 
one parameter. The velocity distribution is given by 


g(v) = A v v ~ a , 


( 4 ) 


To predict the evolution of an electron beam travelling through 
the heliosphere we will rely on the quasilinear approximation 
of wave-particle inter actions (e.g. Drummo nd & Pin es 1962; 
|Vedenov et al. [196 2). This utilises the WKB approximation 
where we are treating waves as quasi-particles interacting res¬ 
onantly (o pe (r ) = kv , where co pe (r ) is the background plasma 
angular frequency (Bian et al.]20T4|. Using /(v, r, t) as the elec¬ 
tron distribution function and W(v, r, t) as the Langmuir wave 
spectral energy density in the one dimension of propagation, we 
can describe the time evolution using 


where a is the velocity spectral index. We have chosen a to be 
7 that corresponds to 3.5 in energy space, and is indicative of 
injected electron beam spectral indices derived from X-ray ob¬ 
servations during flares (e.g. Dennis|1988[|Holman et al.|2011 ). 
The velocity normalisation constant is given by A v = nt, eam (a - 


l)/ v mi~n The parameter nteam is the total density of electrons 
injected at the centre of the injection site between the minimum 
and maximum velocities v min = 2.6vj e = 1.43 x 10 9 cm s -1 and 
\ max - 36 \Te = 2 x 10 10 cm s _1 . We set n\j eam = 2 x 10 6 cm -3 . 

The spatial distribution is given by 



v d 
M(r ) dr 


M(r)f = 


4n 2 e 2 d 

nil 

4nn e e 4 


Wdf 
v <9v 


a f 


mi 


In A— — +S(v,r,f) 

OY V z 


( 1 ) 


dW d(j) L dW dco pe dW _ m» pe ^ df 
dt dk dr dr dk n e dv 

-(n + 7c)W + e 2 (o pe vf In —. (2) 

vr e 

where ojl - (o pe + 3k 2 v Pe /(2oj pe ) is the Langmuir wave dis¬ 
persion relation, m e is the electron mass, e is the electron charge, 
n e (r) is the background plasma density (defined in Section |T.2 .2) 
and In A is the Coulomb logarithm that is assumed constant and 
= 20. M(r) is the cross-sectional area of the magnetic flux tube 
(defined in Section [ 3 ] In this case where M(r) oc r 2 it corre¬ 
sponds tospherically symmetric expansion of the flux tube (e.g. 
|Reid & Kontar|20L?| . Lor the background electron thermal ve¬ 
locity we have v Te = ^fk^^Jn^ e where k B is the Boltzmann con¬ 
stant. The background plasma frequency is defined as (o pe (r ) = 

sj4nn e e 2 /m e . 

Lor the particles, the physical processes we are considering 
in order of terms from left to right in Equation Q are: change in 
time of the distribution function, transport through space in an 
expanding magnetic field, resonant interaction with Langmuir 
waves, collisional damping from the background plas ma, an d a 
source term of electrons S (v, r, t) is defined in Section |2l 2. 1| 

Lor the Langmuir waves, the physical processes we are 
considering in order of terms from left to right in Equation 
[2] are: change in time of the Langmuir wave spectral energy 
density, propagation of waves through spac^J Langmuir wave 
refraction, resonant interaction with the high velocity elec¬ 
trons, Landau damping from the background Maxwellian plasma 
y L = 1 2 ( 0 pe (v/vj e ) 3 exp (-v 2 /2v^J, collisional damping 

from background ions y c - 7in e e 4 \nKI(m 2 e v 2 Te ), and the spon¬ 
taneous generation from the background plasma. 


2.2. Initial conditions 
2.2.1. Electron beam 

We introduce a source of electrons into the simulations, where 
the source function takes the following form: 

S (v, r, t) = g(y)h(r)i{t). (3) 


h(r) = exp [ - ^ ) , 


( 5 ) 


where the characteristic parameter in space is given by the spread 
of the electron beam d. We have chosen d at 10 9 cm which is a 
typical size of a flare acceleration r egion derived from sim ulta- 
neous X-ray and radio observations ( Rei d et al.|201 l|[20l4] ). The 
centre of the injection, r = 0 has an altitude of r in j = 50 Mm, 
consistent with the estimates from a number of flares (Reid 
et al.|[2011, 2014). This altitude corresponds to a density of 
2.14 x 10 y cm 3 an d plasm a frequency of 415 MHz in our den¬ 
sity model (Section 2.2.2| ). Whilst we are injecting quite a large 
number of electrons into the simulation ( nb eam /n e = 10 -3 ), we 
note that the large background density will cause the majority 
of the e lectrons to lose their energy before they can escape the 
corona ( [Reid & Kontar|2013| . 

The temporal profile is given by 


i(t) = A t exp 


(t tin j) 


( 6 ) 


where the characteristic parameter in time is given by the in¬ 
jection time r. We have chosen r = 0.001 s, which is equiv¬ 
alent to an instantaneous injection of electrons in the corona. 
The time normalisation constant is given by A t = l/( y/nr) such 


that i(t)dt = 1. An instantaneous injection was chosen to 
minimise the influence of injection time on the starting height 
that Langmuir waves are induced by the electron beam ( [Reid & 


Kontar |2013| Ratcliffe et al.||2014| ). The constant t in j = 4r is a 


delay time such that four characteristic times can occur before 
the injection maximum. 


2.2.2. Background plasma 

We introduce a population of thermal electrons as a back¬ 
ground plasma. This background Maxwellian population is char¬ 
acterised by a background temperature T e = 2 MK that cor¬ 
responds to a background thermal velocity of Yje - 5.5 x 
10 8 cm s -1 . The choice of 2 MK is related to the higher Landau 
damping that is_predicted from the strahl present in the helio¬ 
sphere (e.g. |Maksimovic et al.|2005| . 

Lor the background electron density n e {r ), we use the den¬ 
sity profile of the Parker model that solves the equations for a 
stationary spherical symmetric solution (|Parker| 195 8[) with nor¬ 
malisation factor found from satellites ( [Mann et al.|1999| . 

r 2 n e (r)v(r) = C = const (7) 


1 Since Langmuir waves do not propagate far before being re¬ 
absorbed, M(r) is not taken into account in Equation (|2j 


y(r ) 2 [n / v(r ) 2 'j 

\T 2 \ \T^ I 


= 41n 



+ 4- 


r 


-3 


( 8 ) 
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where the critical velocity v c is defined such that v c = 
(i k B T sw /pLm p ) 1/2 and the critical radius is defined by r c (v c ) = 
GM. v /2v 2 (both independent on r). T sw is the temperature of the 
solar windrl taken as 1 MK, M s is the mass of the Sun, m p is 
the proton mass and jl is the mean molecular weight. The con¬ 
stant appearing above is fixed by satellite measurements near 
the Earth’s orbit (at r - 1 AU, n = 6.59 cm -3 ) and equates to 
6.3 x 10 34 s -1 . This model is static in time, set at the start of the 
simulations. We justify this because the electron beam is mov¬ 
ing at least two orders of magnitude faster than the solar wind 
velocity. 

We also consider a thermal level of Langmu ir wav es in the 
background p lasma (e.g. |Hannah et al.||20Q9l |2013[ |Reid & 
[Kontar 2013) when wave collisions are weak. This takes the 
form 


W imt (v, r,t = 0) = 


k B T e Upe ^ / V 
An 1 


lnl — 

vr e 


(9) 


where k B is the Boltzmann constant. 


2.3. Langmuir wave energy density 

Usi ng the mode l and the initial conditions described in Section 
|2.1| and Section |2.2| respectively, we have simulated the propa¬ 
gation of an electron beam injected in the corona that can travel 
into interplanetary space. The electron beam resonantly inter¬ 
acts with Langmuir waves forming a beam-plasma structure (e.g. 
|MerNik|1995[|Kontar et al.|1998| ). This structure consists of an 
electron beam inducing Langmuir waves at the front of the beam 
and absorbing Langmuir waves at the back of the beam. The con¬ 
stant generation and absorption of Langmuir waves allows the 
beam-plasma structure to travel large distances up to 1 AU and 
beyond without running out of energy. At any moment in time, 
the majority of energy is contained in the electron population 
and only a small fraction is contained in the Langmuir waves. 

As a proxy for the type III radio bursts we use the energy 
density of Langmuir waves. This is justified because only a small 
fraction of the energy contained in Langmuir waves is actually 
converted to type III radio emission. As such, the bulk of the 
Langmuir wave energy is a good tracer of what times and what 
frequencies the electron beam will be radiating radio waves. 
Moreover, dealing only with the Langmuir waves makes the 
problem computationally tractable for the long simulated dis¬ 
tances of 1 AU.. It has been found that the energy density of 
Langmuir waves is not an exact proxy for the level of radio emis¬ 
sion ( [Ratcliffe & Kontar|2014HRatcliffe et al.|2014[ ) with funda¬ 
mental and harmonic type III radio emission depending differ¬ 
ently upon the spectral characteristics of the induced Langmuir 
waves. However, both branches of emission are dependent upon 
Langmuir waves; the absenece of enhanced Langmuir waves 
will stop both emission types. The radial distance that Langmuir 
waves cease will give a good indication of where the beam in¬ 
stability terminates. No Langmuir waves means no radio waves. 

We calculated the energy density of Langmuir waves E w (r, t) 
(ergs cm -3 ) at every point in time and space by integrating the 
spectral energy density W(v, r, t) over k 

rkmax rVmax qj 

E w (r, t) = I W(k,r,t)dk= I — W(y,r,t)d\. (10) 

Jk min Jv min V 2 

To see the variation of Langmuir wave enhancement over ther¬ 
mal level as a function of distance, we normalised the level of 

2 T sw used in the density model is different from the electron temper¬ 

ature T e that defines v Te 


E w (r , t) by the local thermal level of Langmuir wave energy den¬ 
sity, found from Equation and ( fT0| ) as 


K\r) = 



k B T e te pe 
An 1 v 4 



( 11 ) 


To obtain a stopping frequency for type III bursts, we must quan¬ 
tify at what point in space (and hence what frequency) the elec¬ 
tron beam stops generating high levels of Langmuir waves, re¬ 
quired to emit type III radio bursts. We do this by calculating 
the maximum level of Langmuir wave energy density E™ ax for 
every point of space and time. Again we normalise this by the 
thermal level of Langmuir wave energy density E™ lt . By taking 
a fixed level for E™ ax /E l ” lt we are then able to estimate the low¬ 
est frequency an electron beam was able to induce this level of 
Langmuir wave energy density. We call this the stopping fre¬ 
quency. 


3. Radial expansion of the guiding magnetic flux 
tubes 


When electrons propagate outwards from the Sun, they follow 
the lines of decreasing magnetic field. As the magnetic flux tubes 
expand, so too does the electron beam, causing the beam den¬ 
sity to decrease. An example cartoon showing this behaviour is 
given by Ligure[2](left). We have modelled the expansion via the 
second term in Equation ([!]). We use M(r) to model the cross- 
section of the expanding flux tube that the electrons are travel¬ 
ling along as a function of r and it takes the form 

M(r) = Mojl + Ej . (12) 


where ro is the characteristic length of the expansion of the 
magnetic flux tube. The power-law index defines an expand¬ 
ing flux tube and the radial dependency of the magnetic field 
B oc M(r) _l oc r~P. Mo = 7r(J/2) 2 cm 2 ^ 80 Mm 2 is the cross- 
sectional area of the flux tube at the centre of the acceleration 
site, r = 0. Lor ro = 3.4 x 10 9 cm, an expansion of r 2 gives a 
cone of angle 33 degrees, an angle that has been deduced from 
some in-situ observations of high energy electrons near the Earth 


(e.g.|Lin|1974| |Steinberg et al.|1985| |Krucker et al.|2007[ [Wang 


let al.|2012| ). 

We have not changed the background density profile when 
changing the expansion of the magnetic flux tube. Whislt the 
magnetic field, solar wind speed and background electron den¬ 
sity are all interdependent properties (e.g. solving equations 12- 
14 in |Kontar|2001a| ), it is beyond the scope of the paper to simu¬ 
late how all three change with radial distance together. We focus 
on the effects of the magnetic field. 

As the electron cloud propagates outwards along the expand¬ 
ing magnetic flux tube (Ligure 2]), the electron number density 
decreases with the rate dependent on the flux tube expansion 
rate. Whilst velocity dispersion also causes an increases in the 
derivative df/dv that grows with distance, the decrease in beam 
density from the electrons following the expanding magnetic 
flux tube is greater. These two processes are not the entire pic¬ 
ture but it is clear that the expansion of the magnetic flux tube is 
expected to be a major factor that governs when subsequent type 
III radio emission stops. 


3.1. Simulations 

To illustrate how the expansion of the magnetic flux tube can 
influence the propagation of an electron beam we numerically 
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Fig. 2. A cartoon of an electron cloud (green) travelling upwards from an acceleration region through the corona, following the 
magnetic field. The density of the electron cloud is represented through the transparency. Left: an example of a constant radial 
expansion described by Equation Right: At a critical height r c the magnetic field lines begin to expand at an increased rate 
described by Equation ( f]~3j ), causing the density of the electron cloud to decrease faster. 


solved Equations m with *e parameters given in Section [2] 
except that we change the rate of magnetic flux expansion. For 
this we varied /3 in Equatio n ([12] ) such that /3 = 2, 2.5, 3, 3.5. 

As described in Section |23| we used the Langmuir wave en¬ 
ergy density as a proxy for the type III radio emission. We plot 
the Langmuir wave energy density variation as a function of time 
and frequency (see e.g. Figure [3]), similar to what is normally 
produced for type III bursts using the flux density of radio waves 
as a function of time and frequency. To account for the decreas¬ 
ing Langmuir wave energy as a function of frequency, shown in 
Equation ( fT0| ), we normalised the energy density by the thermal 
level of Langmuir waves given by E™ lt = E w (t = 0) = f Wmdk. 

Figure [3] shows how the energy density of Langmuir waves 
varies as a function of frequency and time for all four values of 
P - 2, 2.5, 3, 3.5. When /3 = 3.5 (top left graph) the energy 
density in Langmuir waves is low, ceasing completely after only 
0.8 mins and not reaching frequencies less than roughly 4 MHz. 
The magnetic flux tube expands very rapidly, causing the elec¬ 
tron beam to decrease in density over a short distance and is con¬ 
sequently no longer able to drive the wave-particle interaction. 
When /3 = 3.0 (top right graph) the magnetic flux tube expands 
less rapidly and Langmuir waves are induced by the electron 
beam over a longer time (2.5 mins) and at lower frequencies (1 
MHz). This trend continues for ft = 2.5 (bottom left graph) and 
P = 2.0 (bottom right graph) with the Langmuir wave energy 
density being significant over the thermal level at 0.1 MHz and 
0.03 MHz respectively. 

For different levels of magnetic flux tube expansion, Figure 
[4] (left) shows the minimum background plasma frequency the 
electron beam was able to excite a c ertain level of Langmuir 
wave energy, as described in Section |2.2.1| We note that the 
background plasma frequency relates to distance using the den¬ 
sity model described in Section 2.2.2 The different lines in 
Figure H illustrate different values of E™ ax IE™ lt , taken between 
10 2 and 10 3 . 

We can see from Figure [4] that slow expansion (smaller p) 
leads to a significant energy of Langmuir waves at frequencies 
below ~ 1 MHz. The points representing different levels of max¬ 
imum Langmuir wave energy density E™ ax /E™ lt are closer to¬ 


gether in minimum frequency for/3 = 2 compared to/3 = 2.5. As 
the electron beam propagates farther from the Sun, the plasma 
density profile becomes flatter, (\doj pe /dr\ decreases) and hence 
the Langmuir wave energy decreases faster with frequency. The 
similarity between the simulations with p - 3.5 and p - 3.0 
is because the electron beam density decreases very quickly in 
both simulations. In comparison to the the cases /3 = 2, 2.5, 
only a small amount of energy is transferred to Langmuir waves 
before the wave generation is quenched. 

When we consider weaker radial expansion (for instance 
P = 2), we observe Langmuir waves being generated at distances 
much farther from the Sun and consequently at much smaller 
frequencies. We considered radial expansion with smaller /3 val¬ 
ues of p - 1.5, 1.0 but they produced so much Langmuir wave 
energy that by 2 AU (the end of our simulation box) there was 
still a significant level of energy density in Langmuir waves, with 
E™ ax / E™ lt > 10 3 . As shown in the next section, part of this result 
is related to not considering the density fluctuatio ns in the back- 
ground plasma. These r esults support the idea (Buttighoffer et al. 
1995; Buttighoffer 1998) that narrow flux tubes are encountered 
when electrons producing type III emission are observed with 
Ulysses at distances greater than 1 AU. 


3.2. Complex Radial Expansion 

In Equation ( [12] ) the magnetic flux tube expands in the corona at 
the same rate as it expands in the heliosphere. Whilst this might 
be more relevant in areas of open field (for instance the coronal 
holes), solar active regions can be far from simplistic. It is possi¬ 
ble that the magnetic field can be focussed in the corona through 
narrow flux tubes. From Section [3] we know this will facilitate a 
higher level of Langmuir wave turbulence. However, it has been 
observed that type III bursts can occur at high frequencies and 
not at lower frequencies. One explanation for this behaviour is 
an increased expansion as electrons propagate out of the corona. 
Figure [2] (right) shows a cartoon of the geometry of the magnetic 
field in such a scenario, where the expansion of the magnetic 
field lines increase at a certain height r c . 
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Fig. 3. Langmuir wave energy density E w IE™ lt in plasma frequency - time plane. The same electron beam parameters are used for 
all panels but with a magnetic flux tube expanding at different rates. The coefficient that models the expansion of the flux tubes is 
P = 3.5 (top left), p = 3.0 (top right), p = 2.5 (bottom left), p = 2.0 (bottom right). Note the different values of the frequency and 
time axis between each graph. 



Expansion Exponent 



4.0 3.5 3.0 2.5 2.0 1.5 

Expansion Exponent 


Fig. 4. The minimum frequency (stopping frequency) where an electron beam is able to induce certain amount of Langmuir wave 
energy density (see color-coded lines), plotted against the different radial expansion exponents, p. The different coloured lines 
represent different levels of Langmuir wave energy (normalised by the thermal level) between 100 and 1000. The right graph 
includes density fluctuations in the background plasma. 


We model the magnetic field proposed in Figure [2] (right) by 
modifying magnetic flux tube cross-sectional area M(r), so that 


M(r) = M 0 


1 + 






for r <r c 
for r > r c . 


(13) 


where r c is a height above the photosphere where the expansion 
of the magnetic flux tube changes. Above this height the expan¬ 
sion is more rapid, with the severity of this increased expansion 


modelled through the constant Z c . Similar to Equation ([12]), the 
constant Mo = 80 Mm is the cross-sectional area of the flux tube 
at the centre of the acceleration site r = 0. 

In Figure [5] we set p = 3 for the magnetic flux tube ex¬ 
pansion, so that magnetic field decreases as r -3 . To simulate 
a reduced expansion in the corona we increased the value of 
r 0 = 10 10 cm. We set r c - 2 x 10 10 cm and r c - 4 x 10 10 cm 
to demonstrate expansion at two different heights above the pho¬ 
tosphere that correspond roughly to 100 MHz and 50 MHz. For 
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Fig. 5. Langmuir wave energy density E w IE™ lt in plasma frequency - time plane. Langmuir waves produced by the electron beam 
with parameters as in Figures [3] [4] but with a magnetic flux tube that undergoes increased expansion at a critical height ( [13] ). Left 
panels: The critical height is r c = 2 x 10 10 cm. The severity of the extra expansion increases from top to bottom. Right panels: The 
critical height is r c = 4 x 10 10 cm. The severity of the extra expansion increases from top to bottom as Z c = 10 1 ,10 3 ,10 5 . 


each height we have varied the constant Z c = 10 1 ,10 3 ,10 5 to 
show different rates of expansion. 

Figure[5]shows the Langmuir wave energy density for the six 
different cases described above. Langmuir wave energy density 
is clearly reduced after the electron beam reaches r c when the ex¬ 
pansion is significant. The two different frequencies of 100 MHz 
corresponding to r c = 2 x 10 10 cm and 40 MHz corresponding to 
r c - 4 x 10 10 cm are evident in Figure[5]when Z c - 10 5 , whereas 
Langmuir waves continue to be induced after the electron beam 
reaches r c when Z c = 10 (Figure [5}. 


Reid & Kontar||2013 ). Langmuir waves are shifted out of res¬ 


onance with the electron beam by the density fluctuations that 
results in a decreased level of Langmuir wave turbulence. 

The power spectrum of density fluctuations near the Earth 
has been observed to obey a Kolmogorov-type power law with 
a spectral in dex of —5/ 3 (e.g. |Celnikier et al.|1983]|1987[ [Chen 


et al.||2013| ). Following |Reid & Kontar | ( |2010| ), the spectrum of 


density fluctuations was modelled with spectral index -5/3 be¬ 
tween the wavelengths of 10 7 and 10 10 cm, so that the new per¬ 
turbed density profile is given by the following equation 


4. Density fluctuations in the background plasma 

The expansion of the guiding magnetic flux tubes is not the only 
factor that determines the stopping distance of a type III burst. 
Density fluctuations in the background plasma are known to de¬ 
crease the level of Langmuir waves induced by an electron beam 


n e (r) = no(r) 


n= 1 


(14) 


(e.g.|Ryutov|1969[|Smith & Sime|1979[|Melrose|1980a[|Kontar| 

|2001c[|Reid & Kontar|2010t|Li et al.|2012t|Ratcliffe et aL|2012[ 


where N = 1000 is the number of perturbations, no(r) is the 
initial unperturbed density (defined in Section 2.2.2| ), A n is the 
wavelength of n-th fluctuation, // = 5/3 is the power-law spectral 
index in the power spectrum, and (j) n is the random phase of the 
individual fluctuations. C(r) is the normalisation that defines the 
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level of density fluctuations 


C(r) = 


2(Aw(r) 2 ) 


5. Role of initial electron beam parameters 

In this section we explore how the different initial electron beam 
parameters can affect the stopping frequency of type III bursts. 
^^ Specifically we look at the number of injected electrons and how 
they are distributed in energy space. 


where the r.m.s. deviation of the density a/(A n{r) 2 ) is taken so 

v 


that near the Earth at r = 1 AU 


= 0 . 1 . 


Reid 


I (An(r=lAU) 2 ) 

(n(r=lAU)) 2 

|& Kontar] ( [2010] ) used numerical simulations of electron beams 
travelling to the Earth to estimate C(r). The best fit was a level 
of fluctuations that decreased as a function of distance using the 
formula 


An _ I (A rc(r) 2 > _ {n 0 (lAU)\' r l(An(r = 1 At/) 2 ) 
n y (n(r)) 2 \ no(r) j y (n(r = 1 AU)) 2 


where 'F = 0.25. This results in a level of fluctuations (denoted 
now for simplicity as An/n) at the Sun that is roughly 1% of 
the level at the Earth, or A n/n = 10 -3 . The value increases as a 
function of distance to An/n = 10 _1 at the Earth. 


4.1. Simulations 

Using the same beam and plasma parameters as Section [3] we 
conducted the simulations with the inclusion of a fluctuating 
background plasma. We varied the expansion of the magnetic 
flux tubes using J3 = 2, 2.5, 3.0, 3.5. The normalised spectral 
energy density of Langmuir waves E™ ax /E™ lt are presented as a 
function of time and plasma frequency in Figure 6] The bursty 
nature of the Langmuir wave energy density is evident, espe¬ 
cially comparing Figure [6] to Figure [3] The bursty nature of the 
Langmuir energy density in Figure [6jis increased for lower fre¬ 
quencies because of the model that we used for density fluctua¬ 
tions. Equation fib] ) has a larger intensity of density fluctuations 
farther away from the Sun. 

Figure [4] (right) shows how the different levels of magnetic 
field expansion affects the minimum frequency of the back¬ 
ground plasma that certain levels of Langmuir wave energy den¬ 
sity are induced. A comparison with the case without fluctua¬ 
tions reveals that in general the inclusion of density fluctuations 
increases this minimum frequency. For instance at = 2.0 an 
increased Langmuir wave energy density was not observed be¬ 
low 0.08 MHz compared to the 0.02 MHz with density fluctua¬ 
tions were not included. Density fluctuations in the background 
plasma are therefore another important parameter in the back¬ 
ground plasma that affects the stopping frequency of type III 
bursts. Increasing the r.m.s. deviation of the density An/n given 
in Equation ( fib] ) increases the suppression of Langmuir waves 
induced by the electron beam density, and vice versa. 

We also point out that for any single value of the points 
in Figure [4] are more closely spaced when you consider density 
fluctuations in the background plasma. Indeed, all points are al¬ 
most on top of each other when= 2.0 or/3 = 2.5. These points 
are due to the spikes in the spatial distribution of Langmuir wave 
energy that reached at least 10 3 times that of the background 
level. The spikes are able to produce higher Langmuir wave en¬ 
ergy density at lower background frequencies. For example a 
level around 10 3 is only observed at 5 MHz without fluctuations 
but at 0.5 MHz with fluctuations. However, the total Langmuir 
wave energy is less when density fluctuations are considered. 


5.1. Injected beam density 

Obviously, the number of electrons we inject into the simulations 
is important when determining the level of Langmuir waves that 
is induced and consequently at what frequency the beam stops 
wave generation. As described in Section [5j the simulation pa¬ 
rameter that dictates the number density of electrons injected 
into the simulation is nt ean 0 Our simulations are set up as de¬ 
scribed in Section [2] with the guiding flux tube characterised us¬ 
ing /3 = 2.5. We have varied the beam density over one order of 
magnitude such that nbeam = 5 x 10 6 cm -3 , 2 x 10 6 cm -3 , 1 x 
10 6 cm -3 , 5 x 10 5 cm -3 . 

Similar to the previous sections, we have plotted the nor¬ 
malised maximum Langmuir wave energy density E™ ax /E™ lt as 
a function of frequency in Figure [ 7 ] As expected, the beams 
with higher initial densities are able to excite higher levels of 
Langmuir wave energy at lower plasma frequencies (or at dis¬ 
tances farther away from the Sun). As the initial number of 
electrons in the beam is reduced, the beam is less and less 
able to generate Langmuir waves farther out. For beam density 
5 x 10 5 cm -3 the beam is not able to excite E™ ax /E™ lt > 10 2,4 
other than when it initially became unstable at frequencies > 
100 MHz, similar to the beams that experience high magnetic 
flux tube expansion in Figure [5] 

5.2. Injected beam energy distribution 

The distribution of beam electrons as a function of energy is very 
important for the production of Langmuir waves. Exactly what 
distribution we inject into the simulation will have a very sig¬ 
nificant effect on the minimum frequency the beam can interact 
with Langmuir waves. The energy distribution is governed by 
the spectral index S defined in Equation 0. We demonstrate 
the importance by using simulations with the same parameters 
as Section |3j keeping the magnetic field expansion static using 
P = 2.5, and varying the velocity spectral index S = 6,7, 8. 

We plot the normalised maximum Langmuir wave energy 
density E™ ax (r)/E™ lt (r) for these three initial spectral indices as 
a function of frequency in Figure [ 7 ] We can see the frequencies 
corresponding to a high peak level of Langmuir wave energy 
density dramatically decreases as we decrease the initial veloc¬ 
ity spectral index. As we have normalised the number density 
of electrons injected, if we lower the spectral index, more high 
energies are injected. There are two reasons this leads to higher 
level of Langmuir wave turbulence: 

- The injection altitude is relatively low in the corona (f pe = 
500 MHz). The background density is 2 x 10 9 cm -3 where 
the collision rate is not negligible. Therefore the low en¬ 
ergy electrons (e.g. < 6 keV) collisionally lose energy be¬ 
fore the beam escapes to a lower density background plasma. 
As demonstrated in [Reid & Kontar| ( |2013| the electrons with 
energy 1.5-6 keV will lose 99% of their energy after only 
0.1 R© of travel. A lower spectral index (harder spectrum) 

3 nb eam is the time integrated beam density injected at the centre of 

the acceleration region. However, given the fast injection characteristic 

time of 10 -3 s, we refer to it simply as the beam density. 
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Fig. 6. Langmuir wave energy density E w IE™ lt in the plasma frequency - time plane. The electron beam initial parameters are as in 
Figures [3] [ 3 ] except that the background plasma includes density fluctuations. Again, the coefficient that models the magnetic flux 
tube expansion is = 3.5 (top left), /3 = 3.0 (top right), = 2.5 (bottom left), = 2.0 (bottom right). Note the different values of 
the frequency and time axis between each graph. 




Beam Density [cm 3 ] Velocity Spectral Index 


Fig. 7. The minimum frequency (stopping frequency) where an electron beam is able to induce certain amount of Langmuir wave 
energy density (see color-coded lines). Left: plotted against beam density. Right: plotted against velocity spectral index. The different 
coloured lines represent different levels of Langmuir wave energy (normalised by the thermal level) between 10 2 and 10 3 . 


means that the beam loses less energy in the low corona and 
retains more energy to transfer to Langmuir waves. 

- Since the growth rate of Langmuir waves at one point in 
space is proportional to v 2 df/dv, the v 2 term means that a 
flatter initial electron spectrum will increase it. 

6. Discussion and conclusions 

To analyse the stopping frequency of type III radio bursts the 
outward propagation of deka-keY electrons have been modelled 


from a flaring acceleration site through the solar corona and 
inner heliosphere. The wave-particle interaction between elec¬ 
trons and Langmuir waves, Coulomb collisions, Langmuir wave 
refraction, and the radial expansion of the guiding magnetic 
flux tubes have been included. We have investigated the low¬ 
est frequency that the beam causes enhanced levels of Langmuir 
waves, and hence radio waves, for different electron beam and 
background plasma properties. 

The radial expansion of the guiding magnetic flux tube sig¬ 
nificantly affects the distances from the Sun that electrons can in- 


9 



























Reid and Kontar: Stopping Frequency of Type III Radio Bursts 


duce high levels of Langmuir wave turbulence (and hence char¬ 
acteristic type III burst frequencies). The electron beam follows 
the expanding magnetic flux tube and consequently decreases in 
density as a function of distance from the Sun. The decrease in 
density eventually quenches the instability that causes Langmuir 
wave generation. This scenario could be the cause of type III 
bursts that are observed only at high frequencies although it is 
not the only mechanism that can produce such type III bursts. 

Background density fluctuations caused an electron beam to 
cease producing Langmuir waves at shorter distances from the 
Sun. However, their effect on stopping frequencies was not as 
significant as varying the rate of magnetic flux tube expansion. 
Due to the clumping of Langmuir waves in space, caused by 
the density fluctuations, we were able to observe high level of 
Langmuir waves at reasonably low frequencies but the overall 
level of Langmuir waves was reduced. 

Varying the properties of the injected electron beam also af¬ 
fects the lowest frequency where Langmuir waves production is 
significant. Decreasing the density of the injected electron beam 
or increasing the spectral index of the injected electron beam 
reduces the distance the electron beam can travel before it stops 
producing high levels of Langmuir wave turbulence. We then ex¬ 
pect that less dense electron beams or beams with high spectral 
indices will have higher type III stopping frequencies. It is pos¬ 
sible that injection of electrons at high altitudes > 0.1 R© would 
be less sensitive to the spectral index of the beam in respect to 
the type III stopping frequency due to a reduced role of Coulomb 
collisions. 

How do these finding relate to some of the hypothesised rea¬ 
sons in |Leblanc et al.[|l995| )? We find that all the different effects 
that govern the number of high energy electrons play some roles 
in the type III stopping frequency: the expansion of the guiding 
magnetic flux tubes, the injected beam density and the injected 
spectral index. We also found that large scale density fluctuations 
in the background electron density play a role, although a high 
level of density fluctuations may still cause irregular radio emis¬ 
sion caused by the localised clumping of Langmuir waves. The 
hypothesis that we have not checked is how the other non-linear 
plasma processes required to convert Langmuir waves into radio 
waves affect the stopping frequency of type III bursts. What we 
have also explained is the prevalence for stronger type III bursts 
to have lower stopping frequencies, a result observed by|Leblanc| 
lefaLl ( [19951 [1995) ; [Dulk et al.| ( [19961 ). The beams with more 
deca-keV electrons produce in general higher levels of Langmuir 
waves, and consequently make brighter type Ills and have lower 
stopping frequencies. The upcoming missions Solar Orbiter and 
Solar Probe Plus should help to observationally study the radial 
plasma properties with in-situ measurements. 

Magnetic flux tubes have been observed to expand super- 
radially (rates faster than r 2 ) in the corona with type III radio 
emission following the magnetic flux tube expansion (e.g. Klein 


|et al.|2008| ). Whilst this expansion can be large, it is very difficult 


to have a large expansion such that B 


.-3.5 


and produce type 


III radio emission. The electron beam decreases in density too 
fast. To demonstrate this we injected an electro n beam with the 
same characteristics described in Section 2.2.1| but with a spec¬ 
tral index of 5 in velocity space. Such a low spectral index is on 
the verge of what is inferred observationally from hard X-rays 
(see e.g. Kasparova et al.|2005| ). Figure [8] shows how the energy 
density of Langmuir waves varies as a function of time and fre¬ 
quency. We only observe a high level of Langmuir wave energy 
at the high frequencies and not at the low frequencies. A flare 
with such a low spectral index would almost certainly produce 
type III radio emission below 0.1 MHz, similar to Figure [I] We 



Fig. 8. Langmuir wave energy density E w IE init in plasma fre¬ 
quency - time plane. The coefficient that models the magnetic 
flux tube expansion is= 3.5 and the velocity spectral index of 
injected electrons is a = 5. 


conclude that we would not observe type III bursts if the solar 
magnetic field ever expands at such a super-radial rate. 

From type III burst simulations ( Ratcliffe et al.||2014| , it 
is evident that the spectrum of Langmuir wave energy density 
does not exactly reflect the spectrum of radio waves that are ob¬ 
served. The requirement for backscattered Langmuir waves to 
interact with forward propagating Langmuir waves in the pro¬ 
duction of harmonic emission means that spectral characteris¬ 
tics of Langmuir waves are important, not just the absolute en¬ 
ergy density. However, the radial expansion of the magnetic flux 
tube reduces the Langmuir wave production. Therefore our con¬ 
clusions on stopping frequency are expected to be mirrored in 
escaping radio emission, although a detailed study including the 
generation of radio emission can explore additional effects of 
wave-wave plasma processes on the type III stopping frequen¬ 
cies. 
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